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Abstract 


A  computer  program  for  calculating  helicopter  torque  coefficients  (Cq)  and 
thrust  coefficients  (Cf)  as  a  function  of  a  vehicle’s  forward  speed  has  been 
developed  in  conjunction  with  the  helicopter  performance  assessment 
project.  The  model  is  based  on  the  energy  principle,  in  which  helicopter 
power  is  broken  into  three  components:  induced  power,  profile  power, 
and  parasite  power.  This  report  documents  the  basic  mathematical  model 
used  in  the  code,  along  with  the  numerical  solution  scheme  used  in 
implementing  the  model.  Results  are  calculated  for  the  UH-60A  Black 
Hawk  helicopter  for  hover  and  different  forward  speed  settings  and 
correlated  with  existing  flight  test  data.  The  effects  of  different  disk 
loading  on  helicopter  power  requirements  are  also  investigated.  The 
present  model  agrees  reasonably  well  with  the  flight  test  data,  providing 
the  author  with  a  certain  confidence  in  the  helicopter  aerodynamic  model 
developed  in  the  present  study. 
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ANALYTICAL  CALCULATIONS  OF  HELICOPTER  TORQUE  COEFFICIENTS  (Cq) 
AND  THRUST  COEFFICIENTS  (Cf)  VALUES  FOR  THE  HELICOPTER 
PERFORMANCE  (HELPE)  MODEL 


1.  INTRODUCTION 

As  a  Fiscal  Year  1998  mission  project,  the  Air  System  Engineering  Analysis  Team  (ASEAT) 
of  the  Survivability/Lethality  Analysis  Directorate  (SLAD)  of  the  U.S.  Army  Research  Laboratory 
(ARL)  was  tasked  to  develop  an  engineering  analysis  model  to  calculate  the  helicopter  steady  state 
performance.  The  model,  which  is  to  be  used  in  support  of  helicopter  ballistic  vulnerability 
assessment  projects,  should  be  capable  of  determining  the  performance  degradation  of  a  helicopter 
when  its  engine  is  damaged  during  combat. 

After  literature  surveys  and  lengthy  discussions  among  team  members  as  well  as  academia 
professionals  (e.g..  University  of  Maryland,  Georgia  Tech  Research  Institute,  etc.)  for  a  possible 
candidate  program,  ARL  decided  to  adopt  the  existing  Helicopter  Performance  (HELPE)  code.[l] 
The  reasons  for  the  selection  of  the  code  are  ease  of  use,  quick  turn-around  time,  and  the 
availability  of  source  code  and  a  documentation  report.  HELPE  requires  many  engineering  data 
to  be  input,  and  most  of  them  are  fixed  parameters  specific  to  the  helicopter  to  be  investigated. 

For  example,  helicopter  gross  weight,  empty  weight,  rotor  diameter,  rotor  rpm,  blade 
aerodynamic  coefficients,  available  engine  powers,  etc.,  are  required  for  the  HELPE  input. 

Other  necessary  input  data  for  the  HELPE  code  is  a  set  of  helicopter  torque  coefficients 
(Cq)  and  thrust  coefficients  (Cf)  as  a  function  of  a  vehicle’s  forward  speed  (ji).  However,  flight 
test  data  specifying  these  values  are  not  readily  available  for  most  helicopters.  To  expand  the 
applicability  of  the  HELPE  model,  a  computer  code  to  numerically  calculate  these  input  values 
needed  to  be  developed. 

In  this  report,  formulation  and  analysis  details  of  the  recently  developed  code  for 
calculating  Cq  and  Cj  of  a  generic  helicopter  system  are  presented,  including  brief  descriptions  of 
basic  input-output  relations  of  the  code.  Numerical  results  are  obtained  for  UH-60  Black  Hawk 
helicopter  using  the  present  code  and  correlated  with  flight  test  data.  [2]  Based  on  results 
obtained,  conclusions  are  drawn,  with  some  recommendations  for  future  works. 


2.  NOMENCLATURE 
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Blade  chord 

Lift  coefficient  at  zero  angle  of  attack 
Lift  curve  slope 
Blade  section  drag  coefficient 
Blade  mean  drag  coefficient 

Blade  section  lift  coefficient 

Blade  section  pitching  moment  coefficient  about  aerodynamic  center 

Fuselage  drag  coefficient 

Fuselage  rolling  moment  coefficient 

Fuselage  lift  coefficient 

Fuselage  pitching  moment  coefficient 

Fuselage  yawing  moment  coefficient 

Helicopter  power  coefficient 

Helicopter  torque  coefficient 

Helicopter  thrust  coefficient 

Helicopter  weight  coefficient 

Viscous  drag  coefficient 

Pressure  drag  coefficients 

Fuselage  drag 

Pitching  moment  coefficients 

Vertical  distance  from  helicopter  center  of  gravity  (CG)  to  hub  center 
Mach  number 

Incidental  Mach  number  normal  to  chord 
Number  of  blades 
Blade  radius,/ 

Rotor  thrust,  lb 

Helicopter  forward  speed,  ft/sec 

Hub  center  position  relative  to  helicopter  CG  in  the  X  and  Y  directions,  respectively 
Distance  between  main  rotor  hub  and  tail  rotor  hub 
Fuselage  side  force 
Tail  rotor  thrust 

Blade  section  angle  of  attack,  radians  (rad) 

Hub  plane  tilt  angle  relative  to  flight  direction 
Longitudinal  shaft  tilt  relative  to  wind  axis,  rad 
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Po  Rotor  coning  angle,  rad 

PlsPlc  Lateral  and  longitudinal  rotor  disk  tilt  angle,  respectively,  rad 

Qlcdjs  Lateral  and  longitudinal  cyclic  trim  input,  respectively,  rad 

0.75  Collective  blade  pitch  at  75%  radius,  rad 

Op p  Helicopter  flight  path  angle  relative  to  the  longitudinal  axis 

Otr  Tail  rotor  collective  control  setting 

jl  Advance  ratio,  V/QR 

a  Rotor  solidity  ratio,  NbCm/nR 

<{>s  Lateral  shaft  tilt,  rad 

p  Air  density,  slug/ft 3 

Q  Rotor  rotational  speed,  rad/sec 

3.  FORMULATION 

3.1  Equivalency  of  Cq  and  Cp 

Let  Q  and  P  be  the  torque  and  power  required  for  a  helicopter  to  sustain  the  level  forward 
flight,  respectively.  Then,  based  on  the  fundamental  physics,  we  have  the  following  relationship 
between  P  and  Q: 

P  =  QO.  (1) 

By  expressing  Q  and  P  in  terms  of  its  coefficients  and  other  helicopter  parameters, 

Q  =  CQpAR(ClR)2  (2) 

P  =  CPpA{OR)3 

Substituting  Equation  (2)  into  Equation  (1),  we  can  obtain  the  following  equivalency  of  Cq  and 
Cp: 


Cq  -  Cp  (3) 

Therefore,  we  can  calculate  the  power  coefficients  instead  of  the  torque  coefficients  for  the 
HELPE  code  input. 

In  the  next  section,  formulations  and  procedure  to  calculate  the  helicopter  power 
coefficient,  Cp,  are  presented. 
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3.2  Helicopter  Power  Coefficient.  Cp 


From  the  well-known  energy  principle  for  rotary  wing  aircraft  [2],  we  can  break  the 
helicopter  power  into  three  fundamental  components:  induced  power  (Pi),  profile  power  (P0), 
and  parasite  power  (Pp). 


P  =  Pi  +  P0  +  Pp  (4) 

In  terms  of  coefficients,  helicopter  power  coefficient,  Cp ,  can  be  expressed  as  follows: 

Cp  =  CPi+CPo  +  CPp  (5) 

in  which  Cpi  ,CpQ ,  and  Cpp  are  the  coefficients  of  induced  power  (Pi),  profile  power  (PQ),  and 
parasite  power  (Pp ),  respectively. 


3.2.1  Induced  Power,  Cpt 

Induced  power  is  the  power  required  to  produce  rotor  thrust.  Using  the  momentum  theory 
[3],  the  induced  power  coefficient  can  be  simplified  as  follows  (for  details  of  derivation,  see 
Reference  3): 

Cp,  =  kX{Ct  (6) 


in  which  Xj  is  the  induced  velocity  coefficient,  Cpis  the  rotor  thrust  coefficient,  and  k  is  an 
empirical  factor  to  cover  nonuniform  flow  and  tip  loss.  Typically,  a  value  of  1.15  is  used  for  k. 
Therefore,  for  given  sets  of  rotor  thrust  (Cp )  and  forward  speed  (ji),  the  only  parameter  needed 
to  calculate  Cpi  is  Xj. 


In  general,  the  rotor-induced  velocity,  Xj,  can  be  calculated  iteratively  using  the  following 
two  equations: 


A,-  =  A  -  /t  tanas 


(7) 


in  which  X  is  the  in-flow  ratio.  In  case  of  hover  (i.e.,  ji  =  0),  Xj  and  X  become  identical: 


(8) 
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Therefore,  the  calculation  of  A,/  requires  another  parameter,  the  longitudinal  tilt  of  the  rotor  shaft, 
as.  This  involves  solutions  of  nonlinear  systems  of  equations  considering  the  vehicle  equilibrium 
condition.  The  formulation  and  solution  of  the  vehicle  equilibrium  equations  are  discussed  in  a 
later  section. 

3.2.2  Profile  Power,  Cp0 

Profile  power  is  a  power  required  to  turn  the  rotor  in  air  (i.e.,  viscous  drag).  Typically,  the 
calculation  of  profile  power  requires  a  knowledge  of  the  blade  section  profile  drag  coefficient, 
preferably  including  its  dependence  on  angle  of  attack  (a)  and  Mach  number  (M).  Often,  it  is 
difficult  to  obtain  a  complete  and  reliable  set  of  even  static,  two-dimensional  airfoil  data. 

Alternatively,  the  rotor  analysis  can  use  a  mean  profile  drag  coefficient  (Cd0)  to  represent 
the  overall  effects  of  the  blade  drag  on  the  rotor.  [4]  The  use  of  a  mean  drag  coefficient  is 
sufficiently  accurate  for  some  purposes,  such  as  preliminary  design,  or  when  detailed 
aerodynamic  characteristics  for  the  blade  section  are  not  available. 

One  of  the  major  setbacks  using  the  mean  drag  coefficient  is  that  rotor  analysis  using  the 
mean  drag  coefficient  often  fails  to  capture  stall  and  compressibility  effects  in  forward  flight. 
Additional  corrections  or  a  more  detailed  analysis  thus  are  required  for  rotors  at  extreme 
operating  conditions. 

In  the  present  calculation,  an  improved  approach  to  calculate  a  mean  profile  drag  coefficient 
using  a  drag  polar  of  the  form  [5]  is  adopted  to  calculate  the  profile  power,  along  with  the  use  of 
Prandtl-Glauert’s  compressibility  correction. 

In  Reference  5,  an  estimation  of  Cd0  is  given  when  the  blade  drag  coefficient  is  known  in  the 
form  of  Cd  =  80  +  SjCC  +  52oc2: 

Q0  =  $o  +  S1CL  +  $2CL 
in  which  CL  is  the  mean  disk  loading: 


c,=^  (10) 

L  O 

Then,  the  coefficient  of  profile  power  is  calculated  as  follows: 
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(11) 


Ce.=?f-{l  +  4.6n2) 


3.2.3  Parasite  Power,  Cpp 

Parasite  power  is  a  power  required  to  overcome  the  drag  of  the  helicopter.  The  estimation 
of  helicopter  parasite  power  is  an  important  aspect  of  performance  calculation  because  it 
establishes  the  propulsive  force  and  power  requirement  at  high  speed.  The  power  required  to 
overcome  the  parasite  drag  of  helicopter  (D),  can  be  expressed  as 

Pp-DV  (12) 


in  which  V  is  the  speed  of  helicopter. 


By  rewriting  Pp  in  terms  of  its  coefficients  and  other  helicopter  parameters,  similar  to 
Equation  (2),  we  have 

C  - 

Pp  pA(QR)3 


(13) 


Further,  the  drag  force  acting  on  a  helicopter  can  be  expressed  as 


(14) 


in  which /  is  the  parasite  drag  area  of  helicopter  including  rotor  hub.  Except  for  compressibility 
or  Reynolds  number  effects,/  is  generally  independent  of  speed. 

By  substituting  Equation  (14)  into  Equation  (13),  we  obtain  the  following  equation  for  the 
parasite  power  coefficient: 


1  _zl_  L=ILus 

2  (OR)3  A  2  A  ^ 


(15) 


in  which  f/A  is  the  equivalent  flat  plate  area  correlated  with  rotor  area  (A).  It  is  found  that  f/A  can 
be  as  high  as  0.025  for  old  helicopters,  about  0.01  to  0.015  for  helicopters  in  current  production, 
and  0.004  to  0.008  for  a  clean  helicopter  design. 

In  this  section,  it  was  shown  that  power  coefficient,  Cp,  can  be  calculated  from  Equations 
(6),  (1 1),  and  (15).  However,  as  mentioned  earlier,  the  solution  of  Equation  (6)  requires  a 
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knowledge  of  rotor  shaft  tilt  angle,  as,  for  a  given  set  of  thrust  coefficient  and  forward  speed,  and 
thus  a  vehicle  trim  solution.  In  the  next  section,  a  formulation  and  solution  procedure  for  vehicle 
trim  is  presented. 

3.3  Vehicle  Trim  Solution 

Propulsive  trim,  which  simulates  an  aircraft  free  flight  condition,  is  used  to  calculate  the 
initial  rotor  control  settings  as  well  as  the  vehicle  orientation  (i.e.,  a5  and  The  solution  is 
determined  from  the  overall  equilibrium  equations:  three  force  (vertical  Z,  longitudinal  X,  and 
lateral  Y)  and  three  moment  (pitch,  roll,  and  yaw)  equations.  These  are 

Fj  =  Fxo  +  DFcosaHP  -  Tsinas  =  0 
F 2  =  Fyo  +  Yfcos([)s  -  Ytr  +  Tsirt(j)s  =  0 
F3  =  Fgg  -  Tcosascos<ps  -  DFcos<j>ssinaHp  +  YFsin((>s  =0 

F4  =  Mxo  +  MXp  +  T{hcosassin<ps  - ycgcos<t>s )  +  YF{hcosascos<ps  +  ycgsin(f>s )  =  0 
F5  =  Myo  +  Myp  +  T[hsinas  -  xcgcosas  ^  +  DF  [~hcosaHP  +  xcgsinaHP )  =  0 

F6  =  Ml+MZp-XtYtr=0 

in  which  Fj,  F2,  and  F3  are,  respectively,  the  force  equilibrium  equations  in  the  X,  Y,  and  Z 
directions,  and  F4,  F5,  and  F5  are  the  rolling,  pitching,  and  yawing  moment  equilibrium  equations, 
respectively.  Also,  Dy  is  the  fuselage  drag;  Yp  is  the  fuselage  side  force;  Ytr  is  the  tail  rotor 
thrust;  T  is  the  main  rotor  thrust;  xcg  and  yCg  and  h  are,  respectively,  the  relative  location  of  the 
rotor  hub  center  with  respect  to  the  vehicle  CG  in  the  X,  Y,  and  Z  directions;  Xt  is  the 
nondimensional  length  (divided  by  rotor  radius)  between  the  main  rotor  hub  and  the  tail  rotor 
hub;  and  as  and  <ps  are  the  longitudinal  and  lateral  shaft  tilts,  respectively.  Furthermore, 

&s  =  V-HP  -Qfp  (1^) 

in  which  oc hp  is  the  hub  plane  tilt  relative  to  the  flight  direction,  and  dFp  is  the  flight  path  angle 
relative  to  the  longitudinal  axis. 

For  a  specified  thrust  coefficient  Cj  and  a  forward  speed  (ji),  the  unknown  quantities  to  be 
determined  from  the  vehicle  equilibrium  equations  are 
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(18) 


u  ^ocs,(f>s,dj^,djC,6]s,Ytr 

These  values  are  recalculated  iteratively  using  Newton’s  method. 

ui+i  =  Uj  +  Auj  (19) 

More  details  about  trim  analysis  and  results  of  several  case  studies  are  available  in  Reference  6. 


4.  RESULTS  AND  DISCUSSION 

Numerical  results  are  calculated  for  the  UH-60A  Black  Hawk  helicopter  for  level  flight 
conditions.  Some  of  important  physical  and  aerodynamic  characteristics  of  this  helicopter  rotor 
are  given  in  Table  1. 

Table  1.  UH-60A  Black  Hawk  Helicopter  Characteristics 


Aircraft  gross  weight 
Number  of  blades,  N& 
Radius,  R 
Blade  chord,  c 
Solidity,  c 

I 

Lock  number,  y 
Blade  airfoils 
Rotational  speed,  Q 
Flap  frequency,  vp 
Nominal  lift  curve  slope,  a 


16,260  lb 
4 

26.83  ft 
1.75  ft 
0.083 
8. 

SC  1095,  SC1095-R8  (tip) 
27  rad/sec 
1.05/rev 
5.73/rad 


Calculated  power  coefficients  for  UH-60A  during  different  rotor  thrust  conditions  are 
correlated  with  flight  data  as  shown  in  Figures  1  through  8. 

Figure  1  shows  the  correlation  of  helicopter  power  coefficients  with  flight  test  data  for  a 
thrust  coefficient  (Cj )  of  0.0066.  The  variation  with  speed  of  the  total  power  coefficient  (C/>  = 
Cpi  +  Cp0  +  Cpp)  and  its  three  individual  components  (i.e.,  Cpi  ,CpQ ,  and  Cpp)  are  plotted  along 

with  the  flight  data.  From  this  plot,  it  is  observed  that  the  induced  power  is  the  largest 
component  in  hover,  but  it  decreases  with  speed.  The  profile  power  exhibits  a  slight  increase 
with  speed.  The  parasite  power  is  negligible  at  low  speeds  but  increases  proportionally  to  /J?  to 
dominate  at  the  high  speed  region.  Therefore,  the  total  power  coefficient  is  high  at  hover  (jj.  =  0), 
has  a  minimum  value  in  the  low  speed  region,  and  then  increases  again  at  high  speed  because  of 
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the  parasite  power.  This  trend  is  quite  agreeable  with  References  2  through  4.  The  results 
obtained  from  the  present  analysis  agree  very  well  with  flight  test,  as  shown  in  Figure  1 . 


UH-60A  HELICOPTER  POWER  COEFFICIENT  VS.  FORWARD  SPEED,  CT=  0.0066 


Figure  1.  Variations  of  Helicopter  Power  Coefficient  (Cp)  Versus  Forward  Speed  (Cr  —  0.0066s). 


Figures  2  and  6  show  the  similar  correlations  but  with  different  rotor  thrust  settings.  Figure 
7  shows  the  correlation  of  power  coefficient  at  very  high  thrust  condition  (Cj  =  0.0102).  In  low 
speed  regions,  the  analytical  results  correlate  well  with  flight  data.  However,  the  present  model 
under-predicts  the  power  coefficients  at  the  high  speed  region  (high  speed/high  thrust).  This  is 
most  likely  because  of  the  shock-boundary  layer  interactions,  which  are  not  accurately  modeled 
in  the  present  study.  In  the  future,  the  use  of  full  Navier-Stokes  computational  fluid  dynamics 
(CFD)  code,  along  with  the  wind  tunnel  experiment,  may  improve  predictions  of  rotor  and 
fuselage  drag. 

The  effects  of  disk  loading  on  the  helicopter  power  required  are  shown  in  Figure  8.  As 
expected,  the  power  required  increases  as  the  rotor  thrust  increases.  This  is  mainly  attributable 
to  high  in-flow  velocity  associated  with  high  disk  loading  at  low  speeds  and  because  of 
compressibility/stall  effects  at  high  speed. 


9 


UH-60A  HELICOPTER  POWER  COEFFICIENT  VS.  FORWARD  SPEED,  CT-  0.0072 
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UH-60A  HELICOPTER  POWER  COEFFICIENT  VS.  FORWARD  SPEED,  CT-  0.009 
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Figure  6.  Variations  of  Helicopter  Power  Coefficient  f  Cp }  Versus  Forward  Speed  (Ct  =  0.0096'). 


UH-60A  HELICOPTER  POWER  COEFFICIENT  VS.  FORWARD  SPEED,  CT-  0.01 02 


Figure  7.  Variations  of  Helicopter  Power  Coefficient  (Cjp)  Versus  Forward  Speed  (Ct  =  0.0102'). 
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UH-60A  HELICOPTER  POWER  COEFFICIENT  VS.  FORWARD  SPEED 


Figure  8.  The  Effects  of  Disk  Loading  on  Helicopter  Power  Coefficient  (Cp). 


Overall,  numerical  results  obtained  from  the  present  analysis  agree  reasonably  well  with 
available  flight  test  data.[l] 

Presently,  another  activity  is  under  way  to  further  validate  the  present  model  with  a  different 
helicopter,  an  AH-64D  Longbow  Apache  attack  helicopter.  Initial  results  from  the  AH-64D  study 
also  show  good  correlation  with  flight  data,  providing  the  author  with  a  certain  confidence  in  the 
helicopter  aerodynamic  model  developed  in  the  present  study. 

5.  SUMMARY 

A  computer  program  calculating  helicopter  power  coefficients  and  thrust  coefficients  as  a 
function  of  a  vehicle’s  forward  speed  is  developed  in  support  of  the  helicopter  performance 
assessment  project.  The  model  is  based  on  the  energy  principle,  in  which  helicopter  power  is 
broken  into  three  components:  induced  power,  profile  power,  and  parasite  power.  For  induced 
power  calculation,  rotor  shaft  tilt  angle  obtained  from  the  vehicle  trim  solution  is  used.  Numerical 
results  are  calculated  for  the  UH-60A  Black  Hawk  helicopter  for  hover  and  different  forward  speed 
settings  and  correlated  reasonably  with  existing  flight  test  data. 
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Based  on  results  obtained  in  this  study  and  further  correlation  study,  the  present  computer 
model  can  be  used  to  estimate  sufficiently  accurate  power  coefficients  for  a  generic  helicopter. 

5.1  Recommendation  for  Future  Study 

In  the  present  study,  an  attempt  was  made  to  include  stall  and  compressibility  effects  in 
the  helicopter  power  calculation  by  using  the  polar  form  of  blade  drag  coefficient  along  with 
Prandtl-Glauert’s  compressibility  correction.  Based  on  results  obtained,  however,  additional 
corrections  or  a  more  detailed  analysis  is  required  for  rotors  at  extreme  operating  conditions  (i.e., 
high  speed/high  disk  loading).  A  future  study  investigating  these  phenomena  would  be  of  great 
interest. 

The  computer  code  can  be  used  to  generate  the  trim  parameters,  as,  6  75, 9]c,  9ls,  required 
as  level  flight  conditions  for  other  ARL  vulnerability  analyses. 
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